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We demonstrate the equilibration of isolated macroscopic quantum systems, prepared in non- 
equilibrium mixed states with significant population of many energy levels, and observed by instru- 
ments with a reasonably bound working range compared to the resolution limit. Both properties are 
fulfilled under many, if not all, experimentally realistic conditions. At equilibrium, the predictions 
and limitations of Statistical Mechanics are recovered. 
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Fundamental aspects of equilibrium Statistical Me- 
chanics (ESM) are intensely reconsidered at present in 
the context of (almost) integrable many-body quantum 
systems [H 0, y, 01 , bringing back to our attention that 
very basic issues are still not satisfactorily understood 
i, i, 0, B i, 0, HH, E d- As in every theory, we 
are faced with the three sub-problems to realistically 
model preparation, time evolution, and measurement of 
a given system. It is well know and will be worked out 
in detail below that the question of experimentally re- 
alistic initial conditions and observables is much more 
urgent in the "derivation" of ESM than in most other 
fields 0, H, H, [m • Regarding time evolution, we 
take the widely (yet not unanimously) accepted view- 
point that standard Quantum Mechanics without any 
additional "postulate" or "hypothesis" must do [isl. [iq. 
The two key questions are then: In how far does a non- 
equilibrium seed evolve to a stationary long-time behav- 
ior ("equilibration")? In how far is this steady state in 
agreement with the corresponding ESM ensemble ( "ther- 
malization")? 

Since open systems (interacting and entangled with 
an environment) are not directly tractable by standard 
Quantum Mechanics, the starting point must be a closed 
(autonomous) system (microcanonical framework), in- 
corporating all relevant thermal baths, reservoirs etc. 

, . Accordingly, the system "lives" in some Hilbert 
space Tl and is at any time instant t > in a mixed state 
(including pure states as special case) p{t) = Utp{0)U} 
with propagator Ut := ex.p{—iHt/h}, seed p(0), and 
time-independent Hamiltonian H. Denoting its eigen- 
functions and eigenvalues by \n) and En (n = 0, 1, 2, ...) 
and the matrix elements {m\p(t)\n) by Pmn{t) we thus 
obtain 

p{t)=Y,Pr.n{0)e^^^--^-^'^'-^\m){n\ , (1) 

where the sum runs over all m, n > 0. As usual, ob- 
servables are represented by Hermitean operators A with 
expectation values Ti{p{t)A} and, without loss of gener- 
ality, are assumed not to depend explicitly on time. 

Generically, the ensemble p{t) is not stationary right 
from the beginning, in particular for an initial condition 
p{0) out of equilibrium. But if the right hand side of 
(H)) depends on t initially, it cannot approach for large 
t any time-independent "equilibrium ensemble" what- 



soever. In fact, any mixed state p{t) returns arbitrar- 
ily "near" to its seed p{0) for certain, sufficiently large 
time-points t, and similarly for the expectation values 
Tr{p{t)A}, see Appendix D in Ref. [13] ■ More specifi- 
cally, consider any p{0) with at least one Pmn{0) 7^ and 
uj := [E„ — E^y^jh ^ 0. Chosing 

A - B + St , S := |TO)(n|/p™„(0) (2) 

it follows that Tr{/9(t)A} = 2 cos(a;t). It is thus clearly 
impossible to "derive" (since it is not correct) ESM for 
arbitrary initial conditions and observables. 

Our first basic assumption concerns the quantity [18] 

R ■■= [ E/^'-W]'/' < [ maxp„„(0)]'/' . (3) 

According to ([T]), the Pnn{t) represent the "occupation 
probabilities" of the energy eigenstates and are inde- 
pendent of t. For a system with / degrees of freedom, 
there are roughly 10'^^-'^^ energy levels in every interval 
of IJ beyond the ground state energy Eq [H, [l^. For 
a macroscopic system with / = 0(10^^), the levels are 
thus unimaginably dense on any decent energy scale and 
even the most careful experimentalist will not be able to 
prepare the system such that the resulting ensemble p(0) 
populates only a few energy eigenstates with significant 
probabilities [I^- For example, assume that there are 
exactly 10'-^° energy levels per J. Even if the system 
preparation defines the energy up to an experimental un- 
certainty of lO-^i^'^^J, there still remain N := loO-S io'' 
energy levels which may be occupied with significant 
probabilities. If all of them are populated equally, we 
obtain p„ri(0) = 1/A^ for N indices n and Pnn(O) = 
for all other n, yielding R < I0~''-3 io according to ([3]). 
If not all N levels are populated equally, but rather any 
p„„(0) may assume arbitrary values between zero and 
10^^" •* times the average population 1/N, Eq. ^ still 
yields R < I0~°-26 io _ Returning to the general case, we 
can conclude that even if the system energy is fixed up to 
an extremely small experimental uncertainty and even if 
the energy levels are populated extremely unequally, we 
still expect that R will be extremely small, typically 

R = 10-°(/) . (4) 

Physical reasons for such a broadly spread energy level 
population include: the time-energy uncertainty relation. 
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imperfect (as opposed to ideal) measurements during the 
system preparation < 0), entanglement processes and 
a time dependence of the Hamiltonian, both caused by 
that part of the environment from which the system is 
isolated for t > but not for i < 0. 

Given p(0), let H+ C H he the Hilbert space spanned 
by those basis vectors \n) for which p„„(0) ^ 0, 

W+ :=span{|7i)|p„„(0) >0} . (5) 

Exploiting Cauchy-Schwarz's inequality [20| 

IPrnn Pinm (t)pnn{t) (6) 

it follows that Pnm{t) = whenever p„„(t) = or 
Pmmit) = 0, i.e. "all non-trivial things are expected to 
happen within . 

Our second basic assumption is that A represents an 
experimental device with a finite range of possible out- 
comes of a measurement within 7^+ [2l| , 

A_A := u\'Ayi{iJ}\A\il}) - min(V'|A|V') = a^ax - amtn , (7) 

where the maximization and minimization is over all nor- 
malized vectors {ip) G Ti.+ and where Umax and amin 
are the largest and smallest eigenvalues of the restric- 
tion/projection of A onto H+. Moreover we require that 
this working range Aa of the device A is limited to ex- 
perimentally reasonable values compared to its resolution 
limit SA, for instance < 10^"°^SA. 

In the worst case, 7i+ — Ti. However, in many cases 
the populations Pnn{0) may be safely negligible e.g. be- 
yond some finite upper energy threshold, yielding a finite- 
dimensional 7i+, while TC is typically infinite dimensional. 
Hence, (O will be finite even for operators A with an un- 
bound spectrum on H. In any case, our above specified 
class of admissible observables A clearly includes any re- 
alistic measurement apparatus. 

For the example ^ we can infer from ([2]) , © , 
that Aa > 2/i?3/2. Hence, the oscillations Ti{p{t)A} = 
2 cos{ujt) are below the resolution limit under our two 
basic assumptions. These (or similar) assumptions seem 
thus indispensable for taming the oscillations in ([T]). 

Given any p{t), we define the auxiliary operator (2^ 

Peq ■=^Pnn{0)\n){n\ (8) 

and focus on the mean square deviation 

al [Tv{p{t)A} - Tv{p,,AW , (9) 

where the overbar indicates an average over all times 
t > 0. The two trace-terms in ([9]) can be unified 
into Ti{p(t)A} with p{t) := p{t) — peq- Introducing 
A :^ A — min-H^ (t/^I A|'0) we can infer from ([7]) that 

< {iJj\A\tI>) < Aa for aU normalized \ip) G H+ . (10) 



Since {n\p{t)\n) = it follows that Tr{p{t)} = and that 
the variance ^ can be rewritten as [Tr{/5(t)A}]2. With 
the help of ([T]) and pnm '■= Pnm(O) we finally obtain 

''A = Y^A.kPkjAranPnm e^[E-E,+E^-EMf^ , (H) 

where the sum runs over all j,k,m,n — 0,1,2,... 
with j ^ k and m ^ n. Next we exploit the fact that 
Ej — Ek + Em — En vanishes for generic 16] Hamiltonians 
H only for j = n and k = m, given j ^ k and m ^ n 
[2^ . Since the time averaged exponentials in (fTT|) vanish 
if Ej — Ek + Em — En 7^ we can conclude that 

aj = ^'|i™„P|p™„P < J2 \^mn\^\Pmn? , (12) 

where the first sum runs over all m n and the second 
over all m, n. With ^ and ([S]) we thus obtain 

fyl !i AjnnPnnA-nmPmm = 

'^{m\Apeq\n){n\Apeq\m) . 

The sum over n amounts to an identity operator and 
that over m yields Tr{[Apeq]^}. This trace over the en- 
tire space Ti. can be restricted to Ti.^ without changing 
its value, as follows from ([5]) and ([5]). Again, this trace 
remains unchanged if we now replace A by Ap := PAP, 
where P is the projector onto H+. Next, we evaluate 
this trace with the help of the eigenvectors |x„) of Ap 
(restricted to H+), yielding 

oi < ^{Xm\PeqAp\Xn){Xn\PeqAp\ 

Observing that Ap\xn) = \Xn) {Xn\Ap\xn) (since \xn) is 
eigenvector of Ap) and that {xn\Ap\xn) = (Xn|^|x«) 
(since |x„) € H+ and thus P\xn) = \Xn)) we can ex- 
ploit (jlOp to obtain 

O-j < A\^{x,n\peq\Xn){Xn\peq\Xm) ■ 

The sum over n yields the identity operator (on 'H+) and 
that over m amounts to Tr{pg^}. With ^ and ^ we 
finally arrive at 

al<AlR'. (13) 

Considering Tr{p(f)^} as a random variable with uni- 
formly distributed i > 0, a similar (but simpler) cal- 
culation as before yields for its mean value the result 
Tr{p(i)A} = Tr{p(.qA}. Hence, dH) is its variance and by 
combining ([T5)) with Chebyshev's inequality [l^], we can 
conclude that 

Prob(^|Tr{p(<)A} - Tr{peg^}| > R A^ < R . (14) 

In view of (|4]) it follows that for the overwhelming major- 
ity of times t >Q the difference between Ti{p{t)A] and 
TrjpegA} is way below the instrumental resolution limit 
5A for any experimentally realistic observable, see below 
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(|7l). In other words, the system looks exactly as if it were 
in the mixed state p^q for the overwhelming majority of 
times t>0 [l^l, though the "true" p{t) is actuaUy quite 
different, see above ([2]). This is our first main result. 

Note that ^4]) is still compatible with the recurrence 
property of Tr{p{t)A} mentioned above ^ but implies 
that such excursions from the "apparent equilibrium 
state" peq must be exceedingly rare events. 

Exactly the same "apparent equilibration" towards peq 
emerges if one propagates p(0) backward in time (keep- 
ing the system isolated). Along the entire real i-axis, an 
initial condition p(0) far from equilibrium thus closely re- 
sembles one of the above mentioned rare excursions, just 
that the location of this excursion is on purpose chosen as 
the time-origin. In other words, the quantum mechanical 
time inversion invariance is maintained, but when start- 
ing out of equilibrium, an "apparent time arrow" emerges 
with extremely high fidelity. 

While p4p provides a bound for the relative amount 
of time the system exhibits notable deviations from equi- 
librium, the typical duration of one given excursion, or 
equivalently, the characteristic relaxation time of an out 
of equilibrium initial condition p(0) remains unspecified. 
Since one can easily find examples with arbitrarily large 
or small relaxation times, any further quantification of 
the relaxation process inevitably would require a consid- 
erably more detailed specification of the Hamiltonian H, 
the initial state p(0), and the observable A. 

Considering and estimating quantities like ^ is very 
natural and has a long tradition: Merits and shortcom- 
ings of the early works are reviewed e.g. in Q, most no- 
tably Ludwig's approach In particular, many of them 
@; 01 involve an extra average over initial conditions with 
the effect that any specific non-equilibrium seed must be 
excluded as "potentially untypical" from the general con- 
clusions. Turning to the more recent precursors, Peres' 
approach 8] is comparable to ours up to Eq. but 
then proceeds with the conjecture that the Amn are pseu- 
dorandom matrix elements, statistically independent of 
the pnrm for which there are general arguments Q and 
numerical evidence \25^ (and counter-evidence [l|) but no 
proof. For pure states, Srednicki obtained similar results 
[ll| by exploiting a common conjecture about the semi- 
classical behavior of classically smooth observables A in 
systems with a fully chaotic classical limit. Again, this 
conjecture is based on good arguments [1^ but no proof. 
Moreover, typical classical many-body systems are not 
expected to behave fully chaotic [IJ. Somewhat similar 
conclusion have been reached even earlier by Deutsch [lOj 
via additional hand waving arguments. Finally, rigorous 
results comparable to ([14]) are due to [3, , but only for 
rather special Hamiltonians H and initial conditions. 

According to ^ and the discussion below pij) , expec- 
tation values become practically indistinguishable from 

Tr{pegA} - ^p„„(0)A„„ (15) 

after initial transients have died out ("equilibration"). In 
how far is this in agreement with ESM, predicting (l5| 



"thermalization" , i.e. the appearance of the microcanon- 
ical ensemble p™"^ instead of peq in In case p™"^ 

and Peq yield measurable differences for experimentally 
realistic p{0) and A, the "purely Quantum Mechanical" 
prediction (115[) is commonly considered as "more funda- 
mental" [ll,U13|- other words, our derivation of ESM 
is complete^ provided the latter is valid itself. In the op- 
posite case, there is nothing to derive, but the quantity 
in (fT5|) still governs the (time-) typical behavior. This 
the second main result of our Letter. 

A first well known validity condition for ESM is a 
"sharp energy E" , i.e. all the /9„„(0) with En G / := 
[E, E + AE] sum up to almost unity, AE being small 
but still experimentally realistic 15J. In particular, the 
Pnn are constant for all En & I and zero otherwise p^ . 
and as a second (often tacit) validity condition for ESM, 
the resulting expectation values Tr{p™"^^} are assumed 
to be (practically) independent of the exact choice of AE 
and E. Basically, this means that the details of pnn{0) 
do not matter in il5\) . henceforth called property (P). 
The same conclusion (P) follows from the equivalence of 
the microcanonical and canonical ensembles (for all en- 
ergies E), considered as a self-consistency condition for 
ESM. Clearly, property (P) is tantamount to replacing 
Peq in by p""'^. Our first remark regarding (P) itself, 
is that no experimentalist can control the populations 
Pnn(O) of the unimaginably dense energy levels En, apart 
from the very gross fact that they are "mainly concen- 
trated within /". If the details would matter, not only 
ESM would break down but also reproducing measure- 
ments, in particular in different labs, would be largely 
impossible. Second, one can readily construct observ- 
ables and initial conditions, being experimentally real- 
istic according to our definitions but still violating (P). 
The fact that ESM is known to have an extremely wide 
experimental applicability implies that our so far notion 
of "experimentally realistic" is still too general. The sim- 
plest option seems to require/assume that the expecta- 
tion values Ann — ('^lAIn) hardly vary within the energy 
interval /. This is similar in spirit to classical coarse 
graining, and, in fact, is part of the already mentioned 
common conjecture about the semiclassical behavior of 
fully chaotic classical systems [2(|. Next, even if the Ann 
notably vary, following Peres '8] the immense number of 
relevant summands in (jl5p may - for "typical" A and p(0) 
- lead to a kind of statistical averaging effect and thus 
a largely p(0)-independent final result. All these conjec- 
tures about "truly realistic experimental conditions" be- 
come even more compelling by considering that, "canoni- 
cally" , A only affects a small sub-system, weakly coupled 
to a "big" rest, which can be readily traced out in (1151) . 
with the effect of an extra averaging step [1, 0, [T^.Tl3j. 
Yet, the apparent universality of property (P) and its 
relation to "more basic" system properties like "ergodic- 
ity" and "(non-)inte grab ility" are still not very well un- 
derstood [E 0, il, [11,111] . 

Numerically, the validity and limits of such conjectures 
and of ESM itself have been exemplified e.g. in [l|, ^ [^, 
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[25j . While the details are not yet settled, "equilibration" 
in agreement with was seen in all cases. Also the 
numerical observation that already quite small particle 
numbers often work surprisingly well is explained by (jl4p 
in view of Uj). 



In th e classical case, proving the counterpart of the 
relation Tr{p{t)A} = Ti{peqA} (see above Eq. fli])) is 
tantamount to the notorious ergodicity problem [3| . The 
next step, namely evaluating the classical counterpart of 
([9]) remains as an even more difficult open problem. 
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